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The  basic  physical  process  considered  is  the  appearance 
of  a  pressure  gradient  in  a  slightly  Ionised  gas  under  the  in¬ 
fluence  of  an  electric  field.  On  &  molecular  scale  this  pro¬ 
cess  involves  the  acceleration  of  ions  in  an  electric  field 
and  the  transfer  of  momentum  to  the  main  body  of  the  fluid 
through  collisions. 

A  particular  phenomenon  involving  thie  basic  process  is  the 
pressure  gradient  in  a  point  corona  discharge  and  the  resulting 
"electric  wind".  The  experimental  arrangement  studied  is  a 
point  to  a  downstream  coaxial  cylinder  discharge  in  air. 

The  difficulity  of  applying  boundary  conditions  to  the 
mathematical  description  of  the  process  requires  some  simplify¬ 
ing  assumptions.  A  simple  one-dimensional  model  of  the  physical 
situation  1b  considered,  the  electrical  field  equations  being  ap¬ 
proached  from  the  consideration  of  a  space  charge  limited  current 
The  functional  relations  derived  appear  consistent  with  the  ex¬ 
perimental  results  for  the  limiting  cases  implied  by  the  assump¬ 
tions. 

Velocity  variations  encountered  in  the  induced  flow  of  air 
present  questions  regarding  the  bulk  parameter  which  describes 
the  molecular  interchange  of  momentum,  the  ion  mobility.  While 
this  experiment  was  not  for  the  purpose  of  measuring  ion  mobili¬ 
ties,  it  is  believed  that  in  any  application's  involving  the  basic 
process  in  the  working  fluid  of  the  aerodynamicist  the  tirae,  den¬ 
sity,  temperature  and  humidity  dependence  of  the  ion  mobility  re¬ 
quires  data  which  is  not  now  available* 
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Current  in  amperes 

Current  density  -  current/area. 
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A  rationalized  length-mass- time-charge  system  is  employed 
where  charge  is  in  coulombs,  and,  compatible  with  the  usual  aero¬ 
dynamic  system,  length,  mass,  and  time  have  units  feet,  slugs,  and 
seconds.  It  is  found  convenient  to  retain  the  volt  as  the  unit 
of  electric  potential  with  the  result  that  force  is  given  in  volt- 
coulombe  per  foot.  In  this  system  the  dielectric  constant  has 
the  value 

€  •  f',70  x  10" 12  coulombs 

volt-feet 

The  hydrodynamic  pressure  in  electrical  units  is  in  volt-coulombs. 

Air  pin  lb  ./ft'", 

1  volt-coulombs  •  . 
ft  3 


7376  lb. 
,ftfc 
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W 1th  a  riev  toward  possible  aerodynamlo  applications  of 
electrical  field  forces  in  an  ionised  gas  there  ie  considered 
here  an  effect  due  to  such  foroee  which,  although  small  by 
comparison  to  the  effects  of  dynamic  forces  generally  conoid* 
ered  ie,  nevertheless,  real  and  measurable*  This  effect, 
which  is  well  known  but  of  only  passing  intereet  to  the  phy¬ 
sicist  and  electrical  engineer,  is  the  ’’electric  wind"  or  the 
flow  of  air  induced  by  a  point  corona  discharge. 

The  complete  problem  involves  an  interaction  between  the 
electrical  field  and  the  flow  field  of  the  air.  The  electric 
field  induces  an  air  flow,  and  the  air  flow,  established  ac¬ 
cording  to  additional  boundary  conditions,  in  turn  affects  the 
electric  field. 

To  prevent  the  analysis  of  the  problem  from  becoming  hope¬ 
lessly  complex  so  that  the  "how  much"  as  well  as  the  "why" 
may.be  answered,  a  simplified  one-dimensional  model  of  the  sy¬ 
stem  ie  considered.  Tor  this  model  apr roximations  ore  required 
only  for  the  boundary  conditions  since  the  governing  field  e- 
ountiona  arc.  quite  simple, 

Initally  considered  is  the  pressure  resulting  from  a  given 
electric  field  and  charge  distribution  -  the  charge  distribu¬ 
tion  being  determined  by  the  properties  of  the  ions,  the  prop¬ 
er'  ies  of  the  air  ( mainly  velocity  and  density),  nnd  by  the 
electric  field.  Secondly,  the  effect  of  the  airstreom  on  the 
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electric  field  le  estimated  from  the  measurable  charac'erls 
tics  of  the  discharge  varying  the  air  velocity*  Then,  for 
given  boundary  conditione  on  the  air  flow  the  induced  airr 
stream  Velocity  may  be  determined  as  a  function  of  the  elec 
trie  current* 
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ii.  £AsssyaE  pua  to  iok  koiios 

Momentum  lacchange  and  Ion  Mobility 

Of  interest  in  this  study  is  the  behavior  of  a  slightly 
ionized  gas  under  the  influence  of  an  applied  electric  field* 
"Slightly  ionized”  refers  to  the  proportion  of  molecules  ion¬ 
ized  and  does  not  imply  that  the  charge  density  or  the  elec¬ 
trical  forces  on  the  charges  need  be  small.  The  condition 
that  the  gas  ia  slightly  ionized  holds  generally  for  a  joint 
corona  discharge  except  for  a  small  localized  region  of  break¬ 
down  ne'-r  the  point.  In  addition,  outside  of  this  localized 
region,  the  ions  are  unipolar,  and  we  may  consider  the  motion 
of  ions  of  only  one  sign. 

On  a  molecular  scale  the  ions  are  accelerated  by  the  e- 
lectric  field  and  interchange  momentum  with  the  neutral  gas 
molecules  b;.  collision. 

On  a  macroscopic  scale  the  bulk  property  which  charac¬ 
terises  this  molecular  interchange  is  the  ion  mobility.  If 
uj  ia  the  average  ion  drift  velocity  and  E  the  magnitude  of  the 
electric  field,  then  the  ion  mobility,  k,  is  defined  by 
Ui  a  kid 

To  show  that  k  io  actually  a  parameter  which  relates  the 
electric  current  to  the  momentum  exchange  with  the  gas,  con¬ 
sider  ion  motion  through  a  stationary  gar.  Multiply  both  sides 
of  the  defining  equation  for  k  by  the  charge  density,  CT  , 
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cr  uj  •  ks<r 


or 
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flow,  <T  ui  is  just  the  current  density,  j,  while  «TE  la  a  force 
per  unit  volume  or  better  a  pressure  gradient  exerted  on  the 
fluid  '>r  conductor. 

The  ion  mobility  ia  readily  related  to  the  conductivity 
of  a  continuous  medium.  If  the  conductivity,  c,  is  defined  by 

t 

the  relation 

j  *  cTC  where  j  ia  the  current  density 

then  c  «  cr  k 

Tlic  conductivity  13  more  familiar  as  describing  the  macrosopic 
properties  of  electron  motion.  In  this  case  the  momentum  trans¬ 
fer  to  a  continuous  medium  is  extremely  small,  k  being  of  the 
order  of  10^  times  the  value  for  that  of  negative  ions  in  air. 

The  functional  relation  of  k  with  the  molecular  properties 
of  the  ions  and  neutral  molecules  may  be  derived  using  the  kin¬ 
etic  theory  of  gases  af J er  some  suitable  choice  of  a  molecular 

model.  Different  molecular  models  result  in  differences  in  the 

(6) 

calculated  temperature  dependence  of  k.  Phillips  measured  the 
mobility  of  positive  p.nd  negative  ions  in  air  at  constant  pres¬ 
sure  for  temperatures  between  94°  and  411°  K»  These  results  re¬ 
duced  to  constant  density  (using  the  perfect  gas  equation)  show 
the  mobility  independent  of  temperature  above  P50°  K.  In  agree¬ 
ment  both  experimentally  and  ns  derived  for  all  molecular  models 
is  the  density  dependence  of  k.  '/’he  mobility  is  inversely  pro-. 
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portional  to  the  density  of  the  gae.  | 

The  experimental  values  of  the  ion  mobility  are  not  too 
well  defined  except  for  well  controlled  experiments  with  pure 
gases  where  there  is  some  certainty  about  the  nature  of  the 
ions.  There  is  especially  some  uncertainty  about  the  vnlue 
of  k  for  the  following  experiments  which  were  performed  using 
room  air. 

In  air  the  nature  of  the  ion  ie  difficult  to  define. 

This  i 8  especially  true  for  the  posit ive  ion  for  which  there 

is  a  characteristic  time  at  which  ions  of  a  different  mobility 

appear.  Experiments  indicate  that  for  small  time  the  mobility 

of  the  positive  ion  in  approximately  the  same  an  that  for  the 

stable  negative  ion  -  the  mobility  decr<. using  with  time.  This 

"aging”  of  the  positive  ion  is  generally  accepted  to  result 

from  combination  of  the  ions  with  other  neutral  molecules. 

To  give  an  order  of  magnitude  of  the  characteristic  time  of  the 

(?) 

positive  ion  Loeb  states  that  "no  more  than  two  mobilities  and 
usually  only  one  mobility  appear  for  ionn  in  a  pure  gas  at  any 
age  below  10-?  seconds." 

The  primary  negative  ion  in  air  in  a  corona  discharge  is 
formed  by  electron  attachment  to  oxygen  molecules.  Thus,  in 
a  negative  point  corona  discharge  there  1b  also  a  characteristic 
time  for  the  formation  of  the  negative  ions.  This  time  dependence 
>f  the  ion  mobility  will  later  be  considered  as  a  possible  ex¬ 
planation  of  variations  in  the  induced  velocity  of  the  air  in  the 
corona  discharge. 

for  certain  limiting  conditions  the  ion  mobility  raay  be 
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considered  constant.  As  an  example*  if  the  electrode  geometry 
and/or  the  airstream  velocity  are  such  that  the  time  of  passage 
of  a  positive  ion  is  less  than  the  characteristic  time  for  which 
ion-molecule  combinations  are  probable  then  the  mobility  may  rea¬ 
sonably  be  assumed  constant. 

t'ith  regard  to  the  density  dependence  of  k,  in  the-  experi¬ 
ments  of  this  study  the  airstrenm  velocities  range  only  to  the 
order  of  fOCft/stc.  allowing  the  approximation  of  an  incompress¬ 
ible  flow  of  air.  It  will  further  be  seen  that  the  electric 
field  forces  are  an  order  of  magnitude  smaller  than  the  maximum 
dynamic  forces  of  the  airstream  so  that  density  changes  due  to 
pressure  may  be  neglected.  Thus,  throughout  the  region  of  the 
discharge  the  ion  mobility  will  be  assumed  independent  of  den- 
b  i  tv  . 

Approximate  values  for  the  "nged"  positive  ion,  k"*’,  and 
for  "new"  positive  ions  and  negative  ions,  k”  ,  in  dry  air  at 
standard  sea  level  density  are 

k +  *  1.5  x  10“ 3  ft/sec 

volts/ft,  (1) 

Refj  Thomson  p.irs 

k~  s  f.S  x  10-?  ft/sec 

voltr./ft . 


pressure  Gradients  In  an  Ionised  gas 

Consider  the  motion  of  a  uniform  airstream  with  velocity, 
u,  in  the  direction  of  a  high  voltage  electrode  point*  Outside 
of  the  small  region  of  break-down  near  a  positive  point  there 
is  a  density  of  positive  charge,  O',  which  moves  under  the  ac¬ 
tion  of  the  airstream  and  the  electric  field,  E,» 

For  this  one-dimensional  incompressible  flow  of  air  the  mo¬ 
mentum  equation  of  continuum  mechanics  reduces  to 
Vp  *  OTP 

i.e.,  the  presoure  gradient  is  proportional  to  the  charge  den¬ 
sity  and  the  resultant  electric  field. 

The  conservation  of  charge  is  expressed  by  the  equation 

T  s  <T  (u+kh) 

7'hcre  j  is  the  current  density  (current/area)  and  k  is  the  ion 
mobility*  ':’hc  sum  (TT  +  kE)  represents  the  resultant  ion  veloc¬ 
ity  —  the  ion  drift  velocity,  kE,  now  being  measured  relative 
to  the  moving  airstream. 

Por  the  stationary  case,  u»0,  the  problem  is  much  simplified 

since 


the  current  now  having  the  same  direction  as  S* 
Thus 


V  p  -  111, 

iTFI 


=  1 

k 


since  as  before  E  and  J  are  collinear 
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This  rather  simple  result  shows  that  for  this  stationary 
case  the  pressure  gradient  may  he  determined  without  any  ex¬ 
plicit  knowledge  of  the  field  strength  or  charge  density  which 

may  he  difficult  to  determine  for  all  hut  the  simplest  of 

(4) 

geometries.  Chattock  used  this  result  to  measure  the  mobil¬ 
ities  of  ions  produced  in  a  point-to.-ring  corona  discharge. 
Enclosing  the  electrode  pair  in  a  glass  cylinder  and  making 
the  distance  between  the  electrodes  large  compared  to  the  di¬ 
ameter  of  the  ring  anode,  the  current  was  assumed  one-dimen¬ 
sional  giving 


and 


iA 

Ak 


where  i  is  the  total  measured  current, A p  the  pressure  dif¬ 
ference  along  the  tube,  A  the  cross-sectional  area  of  the  tube, 
and  &  the  axial  distance  between  the  electrodes,  filth  this 
arrangement  k  was  determined  by  measuring  the  pressure  differ¬ 
ence  and  the  current,  giving  the  "right  answer"  after  some 
manipulation  with  the  geometrical  parameters. 

For  the  case  of  an  air  flow  superimposed  on  the  discharge 
the  problem  is  much,  more  complicated  even  with  the  one-dimen¬ 
sional  assum  tion  which  will  be  maintained.  The  one  dimensional 
immentuai  equation  is 

sffp 

dx 
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where  now 
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giving  d£  a  j  S _ 

dx  u+kE 

In  an  integration  over  x  the  current  density  might  leg¬ 
itimately  be  assumed  constant  if  the  electrode  spacing  is  large 
compared  to  the  radius  of  the  outside  electrode*  However, 
the  field  strength  now  appears  explicitly.  An  attempt  to 
apply,  the  simplifying  assumption  of  a  one-dimensional  field 
leads  to  some  difficulty  with  boundary  conditions  especially 
at  x  *  0,  i.e.,  in  the  neighborhood  of  the  point.  Here  the 
field  is  highly  non-uniform,  a  condition  which  1b  necessary  • 
for  the  occurrence  of  the  corona  discharge. 

The  problem  of  determining  the  pr>  ssure  difference  across 
.J,  the  corona  discharge  from  the  above  equation  resolves  itself  to 

finding  an  approximate  one-dimensional  expression  for  the  elec¬ 
tric  field.  This  requires  a  consideration  of  the  characteris¬ 
tics  of  the  point  corona  discharge  and  how  these  characteristics 
arc  affected  by  an  airs  jam. 
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The  Point  Corona  Discharge 

'The  electrical  discharge  at  a  point  nay  be  described  in 
various  phases.  For  a  positive  point  the  phases  may  be  enu¬ 
merated  ns  follows t 

1«)  Conduction  below  ionization  ty  collision. 

"For  low  applied  voltage  the  current  to  the  point  de- 
.  pends  upon  ions  produced  by  an  external  agent#  In  normal 
room  air  the  magnitude  of  this  current  is  minute  and  dis¬ 
continuous# 

P.)  Conduction  with  Ionization  by  collision. 

For  higher  applied  voltages  the  field  strength  in¬ 
creases  to  the  point  where  electrons  resulting  from  external 
ioniztion  gain  sufficient  energy  in  the  field  to  cause 
further  ionization  by  collision#  For  sufficiently  large 
field  strength  the  "Townsend  avalanche"  effect  results  in 
a  current  pulse,  which  greatly  magnifies  the  effect  of  the 
originrl  electron. 

7> .)  The  Ceiger  Counter  regime. 

For  still  higher  field  strength  a  single  externally 
produced  ion  results  in  an  "avalanche"  in  which  sufficient 
radiation  la  produced  to  cause  further  photoelectric  ion¬ 
ization  in  the  gas.  A  large  current  pulse  builds  up  to  a 
point  where  the  space  charge  distortion  of  the  field  chokes 
the  process  of  ionization  by  collision  and  the  discharge 
censes. 


i  ipn  mpum 
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4, }  Continuous  corona  regime*  f 

| 

For  yet  higher  applied  voltages  the  opace  charge  die-  I 

i 

tort  ion  cannot  reduce  the  field  strength  to  the  point  of  I 

choking  and  there  results  a  self  sustained  discharge*  { 

The  visible  corona  glow  now  appears  about  the  point*  This  j 

discharge  can  be  maintained  in  the  absence  of  another  e-  ] 

lectrode  by  photoelectric  ionization  in  the  gas  and  at 
lower  potentials  by  additional  photoelectric  ionization  at 
another  electrode. 

5. )  Breakdown  stage* 

The  field  strength  may  be  increased  to  the  point  where 
the  breakdown  of  the  gas  extends  across  the  electrode  gap. 

This  may  initially  he  in  the  form  of  ’’pre-breakdown  stream¬ 
ers"  or"brushe8".  Finally,  the  corona  discharge  culminates 
in  an  arc  discharge. 

The  only  phase  of  the  discharge  which  is  considered  here 
is  the  continuous  or  self  sustaining  discharge.  Further,  it  is 
considered  only  from  the  viewpoint  of  continuum  mechanics  with-  . 
out  reference  to  the  "mechanistic”  descriptions  of  electron  av¬ 
alanches,  etc.  From  this  point  of  view  there  are  no  dietinctione 
to  be  made  in  the  description  of  the  electric  field  in  a  posi¬ 
tive  or  negative  discharge.  The  only  difference  considered  is 
that  of  the  bulk  parameter  k,  the  ion  mobility. 
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Point  to  Cylinder  Discharge 

The  geometrical  electrode  arrangement  to  he  considered  is 
that  of  a  point  to  coaxial  cylinder  with  the  oylinder  displaced 
downstream  of  the  point  as  illustrated  in  figure  5.  This  ar¬ 
rangement  presents  itself  as  aerodynamically  “clean*  with  re¬ 
spect  to  the  flow  of  air.  This  is  important  ainoe  in  the  meas¬ 
urements  made  the  electrical  pressure  forces  are  of  the  same  or¬ 
der  of  magnitude  as  the  frictional  forces. 

The  description  of  the  electric  field  in  such  a  discharge 
presents  a  difficult  task  if  an  exact  solution  ie  desired.  Even  j 

for  the  simpler  case  of  no  space  charge,  exact  solutions  are 
readily  obtained  only  for  very  simple  geometrical  electrode  ar- 
rangements,  i.e,  where  the  boundary  conditions  ere  simply  given 
in  an  appropriate  coordinate  system.  With  appreciable  space 
charge,  which  is  the  case  in  the  region  near  the  point  electrode, 
the  boundary  conditions  to  be  applied  in  the  solutions  of  the 
field  eouation  arc  not  only  difficult  to  express  mathematically 
but  are  difficult  to  determine, 

’Reasonably  accurate  solutions  for  the  infinitly  long,  thin, 

wire  to  coaxial  cylinder  corona  discharge  hare  been  developed 

(2) 

and  are  presented  by  Loeb.  A  solution  of  Laplace’s  equation  is 
derived  neglecting  the  space  charge  outside  of  a  cylindrical 
surface  at  which  the  boundary  condition  is  applied.  This  pro¬ 
cedure'  of  neglecting  the  space  charge  and  using  Laplace’s  e- 
'•uation'  an  an  approximation  to  the  field  seems  to  be  a  generally  ■ 
accepted  method  of  approximating  the  field  for  general  electrode 
arrangements. 
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If  this  approximation  is  accepted  for  the  point  corona 
discharge,  then  it  would  be  expected  that  pressures  due  to 
momentum  exchange  would  dominantly  occur  in  the  vicinity  of  the 
point  electrode  where  the  field  forces  ns  given  by  Laplace's  e- 
auation  are  extremely  high,  falling  to  negligibly  small  values 
over  a  diet  nee  small  compared  with  that  of  the  electrode  spac¬ 
ing.  Evidence  that  this  is  not  the  case  is  presented  by  Chat- 
(4) 

tock  in  his  preliminary  investigations  of  the  pressures  on  a 
flat  plate  in  a  point  to  plane  discharge.  Pressure  distributions 
were  obt- ined  at  constant  current  and  electrode  spacing  for  two 
cases.  The  first  case  wr  a  for  a  simple  point  to. plane  geometry. 
The  second  case  involved  an  additional  electrode  surrounding  the 
joint  which  collimated  the  field  except  in  the  neighborhood  of 
the  point.  With  the  field  collimated  the  preseure  distribution 
v/r.s  sharply  peaked  just  opposite  the  point  and  fell  off  rapidly 
with  the  radial  dist  nee  from  the  centerline.  Without  the  colli¬ 
mator  the  pressure  distribution  w  a  relatively  flat  with  appreci¬ 
able  pressure  being  maintained  on  the  plate  to  -he  limit  of  the 
traverse  mechanism.  Chattock's  conclusions  were  that 

"The  pressure  cannot  be  due  either  v/holly  or  in  considerable 
part  to  the  imparting  of  momentum  to  the  gas  at  the  point  itself; 
otherwise  the  pressure  curves  a  and  b  would  be  alike,  since  the 
field  at  a  discharging  point  is  constant  for  a  given  current". 

This  evidence  and  the  apparent  success  obtained  by  Chattook. 
in  using  a  one-dimensional  analysis  in  the  point  to  ring  discharge 
discouraged  the  use  of  an  electric  field  calculated  by  neg¬ 
lecting  space  charge  and  applying  this  result  to  determining  the 
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pressures.  Such  an  analysis  would  be  productive  if  a  suitable 
boundary  condition  could  be  calculated  and  applied  at  some  sur¬ 
face  outside  of  the  immediate  neighborhood  of  the  point  i,e», 
at  some  surface  outside  of  which  the  space  oharge  might  legit¬ 
imate  ly  be  neglected  in  describing  the  field*  Such  an  analysis, 
while  complex  for  a  stationary  fluid, bee -mes  further  complicated 
for  the  present  case  in  which  the  air  stream  and  the  discharge 
are  interdependent.  The  airatream  influences  the  boundary  con¬ 
ditions  of  the  electric  field  and  the  electric  field  in  turn 
influences  the  induced  velocities  due  to  the  pressure . gradient. 

The  simplified  model  to  be  considered  is  that  of  a  onefdim- 
ensional  space  charge  limited  current  in  a  one-dimensional  in¬ 
compressible  flov;  of  air.  The  one-d'imens ional  assumption  i8 
maintained  throughout,  including  the  ion  source.  The  ions  are 
assumed  to  emanate  from  a  planar  sourc--  at  a  potential  equal  to 
that  of  the  starting  otentinl  of  the  continuous  corona  discharge. 
The  potential  drop  or  anode-fall  is  taken  as  the  difference  be¬ 
tween  the-  ar.pli<d  t  otentinl  and  the  'at'  rting  potential  of  the 
discharge.  This  energy  loss  maybe  considered  as  that  required 
for  creation  of  the  ions  including  loanee  by  radiation,  excited 
:olec-;lar  states,  and  chemical  reactions  as  well  vo  the  energy 
manifesting  itself  in  useful  ionization. 
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The  Electric  Field  in  One  Dimension 

If  D  is  the  electric  excitation  in  unite  charge  per  unit 
area  and  o’  le  the  charge  density  in  units  charge  per  unit  vol¬ 
ume,  then 

V' D  *  S’ 

TTsing  the  electrostatic  charge-force  relation  given  by 

where  t  is  the  dielectric  constant!  the  differential  equation 
for  the  electric  field  vector,  E,  is 


For  the  one-dimensional  cnee 

dE  -  SL 
dx  “  fc 

hsing  the  equation  expresning  the  conservation  of  charge, 

this  becomes 

dK _ i_  .1 

dx  “  t  u  -rkE 

or  . 

1  d(u  +  kI3)f=i.  dx 

w- 

Integrating  between  the  electrodes  from  x  ■  0  which 
is  taken  to.  be  the  edg<..  of  the  small  region  of  break-down 
about  the  point  anode  and  letting  J30  be  the  field  strength 
x  *  0 

(u+k;-)?  -  .(u  +  k!U0)r  «  r  kjx 
”  *  1  F  2  kjx  +  (u-hkEQ)2|  *  -  u 

s  [r  J  F 
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This  expression  stay  be  further  simplified  if  the  one* 
dimensional  approximation  Is  maintained  by  assuming  that  the 

U># 

Ion  source  la  a  planar  ionised  layer}  then  Thompson  shows 
that  the  field  strength  at  the  edge  of  euch  a  planar  source 
my  be  neglected.  !?ith  this  assumption  the  expression  for 
the  electric  field  becomes 


TTnving  an  expression  for  the  field  strength,  the  equation 
for  the  pressure  gradient  in  terms  of  J  and  u  may  be  written 
down  immediately.  Trowever,  from  the  pressure  gradient  it  would 
be  difficult  to  experimentally  check  the  validity  of  the  one* 
dimensional  assumption  ove any  rcanonable  range  of  air  stream 
velocities  since  electric  field  forces  become  very  small  in  com¬ 
parison  to  the  dynamic  preonure  forces  and  the  viscous  forcce 
at  the  boundary  of  the  flow. 

The  range  of  velocities  and  currents  for  which  this  t>« 
uation  may  be;  valid  will  be  estimated  by  e nrry'ng  the  one*dim* 
eneional  analysis  further  to  obtain  a  relation  for  the  current 
in  terms  of  the  velocity  nnd  applied  potential,  Theoe  relat¬ 
ione  are  r<a  ily  checked  experimentally. 


#  p.  £06 
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One-Dimensional  Voltage-Current  Relation 

As  an  electrostatic  problem  the  field  ie  irrotational,  i*®* 
S  *  -  djfc  where  $  is  the  electric  potential. 

Thus, 

s  •  s  dx 

Integrating  the  expression  for  E  from  x«0  to  x»A  and  letting 
A4>  a  $<,  “  4>i, 

A  $>  *  J  Edx  ■  ^  j*  jjajkx  +  u2J  ^  -  u"^  dx 


•  2  1 f  fajkjt  u?l  5"  -  u3*l-iu  where  a  •  2 

3  ajks  L  -*  if  t 

To  get  an  approximate,  explicit  expression  for  the  current 
the  magnitudes  of  the  various  terms  must  be  evaluated,  especially 
the  constants  of  integration  The  current  as  a  function  of 

applied  potential  is  zero  up  to  a  potential  V  *  V*  which  ie  de¬ 
fined  as  the  starting  potential  of  the  self  sustaining  oorona 
discharge.  The  current  is  primarily  a  function  of  (V  -  V8), 
the  difference  between  the  applied  potential  difference  and  the 
starting  potential. 

In  an  attempt  to  avoid  the  details  of  the  non-uniform  field 
assume  that  A$*  V  -  V8,  i*e,,  that  the  anode  fall  of  potential 
is  V  -  V8.  Considering  first  the  simpler  case  for  u  ■  0 


A  $  .  V 


Vs  •  ||  [ajki.]  * 


Before  checking  the  validity  of  this  equation  in  the  iight  of  ex¬ 
perimental  results  the  order  of  magnitude  of  the  quantity  ajkjfc  , 
which  has  the  dimensions  pf  a  velocity  squared,  will  be 
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t 

estimated.  This  velocity  is  in  the  one-dimensional  case  for 
u  m  0  the  maximum  ion  drift  velocity  attained  at  the  down¬ 
stream  electrode.  Corresponding  to  experimental  values »  leti 

Si  t  axial  distance  between  electrodes  •  6*25  x  10"2  ft. 
A  m  cross  sectional  area  ■  2.08  x  10*“  ft. 

In  addition*  let 

a  m  7.41  x  1011  volt-ft  -w  dielectrio  constant  in  vacuo 
coulomb 

J  ■  i  x  10" 6  where  i  is  the  measured  current  in  micro- 

I 

amps. 

arid  assume 

k  ■  1.6  x  10“ 3  ft/sec 

■*).  volte/ ft 

then, 

aJkJL  *(190  ft  //«.amp. 

Bee 

and  the  current  in  microamps,  is 
i  -  4.18  x  10"8  (V-V8)2 

The  current  -  voltage  relation  of  the  point  corona  is 
usually  described  by  an  equation  of  the  form 
i.CV  (v-v8) 

This  agrees  with  experimental  results  over  a  large  range  of 
voltages,  TTowever,  for  large  currents  the  data  is  better  de¬ 
scribed  by  a  relation  of  the  form 
i  s  C  (V-V8)2 

--  7. or  large  currents  where  the  approximation 

ajk£  **  (k~0) *  is  valid  the  experimental  results  cor¬ 
responding  to  the  above  equation  gives 
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1  ■  5.S  x  10"8  (V-Va)2 

The  agreement  with  the  approximate  equation  is  quite 
satisfactory  since  the  concern  here  is  primarily  to  establish 
some  justification  for  approximate  relations  for  the  pressure 
difference  across  the  discharge.  Better  agreement  in  the 
constant  of  proportionality  might  be  obtained  by  some  manipu¬ 
lation  of  the  geometrical  parameters.  For  instance  the  dis¬ 
tance,  l,  may  be  considered  only  as  a  characteristic  length, 
rsing  the  distance  between  the  electrodes  for  this  length  is 
correct  only  as  a  orientating  value. 
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Current  as  a  "function  of  Voltage  and  Airstream  Velocity 

TTnder  the  influence  of  a  moving  airstream  the  current  ie 

given  implicitly  in  the  one-dimenBional  case  by 

V  -  V8  ■  a  1  T  P ajkjl  +  u2l  «T  -  u3  1  -  u  JL 

3  ajk*  u  J  T~ 


?or  aikl 


Vs  «  2  1  1  |*(ajkj£  )£  +  3  u2  (ajk£  )^+“l  -us[-u£ 

7  aJF"  \  L  r  J  J  r 


j-  2  J  (ajki  £  u2  ,  -  2  £  u3  -  u  1 

"u  +  Tc  Tajrrr*  sssiEr  r 

Introducing  the  velocity  ratio  parameter  Ur  ■  u 

(aJkJ&  )t 

v-vs«||  (ajk£  -  u|[l-ur+|ur2] 
"Dividing  by  the  current  density  at  zero  velocity,  J0,  and 


rearranging 


•  i  +  u_jl 

IQ 


Volk 


(1  -  ur  g  ur2) 
3 


This  relation  expresses  the  characteristics  of  the  dis¬ 
charge  in  an  airstream.  The  major  physical  argument  leading  to 
its  derivation  is  that  the  ions  enter  the  electric  field  with 
the  velocity  u  and  the  electrostatic  potential  energy  (V-V8), 
The  equation  is  still  implicit  in  i  since  ur  is  a  function 

of  i.  A  simple  explicit  expression  is  obtained  for  the  case 
„  « 1 ,  viz., 


[il*  s  1  a  u  JL 
lioj  (V-Vs)k 


(U  T<<1  ) 


This  relation  expresses  the  experimental  results  for  a 
limited  range  of  voltage  and  velocity  implied  in  the  derivation, 
"^or  high  applied  voltage,  i,e,,  with  increasing  space  charge, 
tMs  relation  appears  to  apply,  ao  can  b.e  seen  from  figure  f. 


h 
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nihe  deviations  from  the  linear  relation  of 


occur 


with  decreasing  applied  voltage. 

In  attempting  to  evaluate  the  geometrical  factor,  Ji  ,  using 
the  supposedly  more  accurate  expression  involving  the  expansion 
in  terms  of  the  velocity  ratio,  ur,  it  was  found  that  the  ap¬ 
proximate  expression  for  the  current  was  a  better  functional  re¬ 
lation  in  u.  ’’’his  would  further  limit  the  use  of  the  assumed 
expression  for  the  field  strength  to  the  condition  that  ur«  l. 

The  data  of  figure  f  represents  airatream  velocities  only 
as  low  as  75  ft/sec.  figure  3  indicates  that  the  data  may  be 
line-rly  extrapolated  to  lower  velocities. 

7rom  the  above  it  is  expected  that  the  equation  for  the 
assumed  electric  field  may  be  a- plied  to  evaluating  the  pressure 
under  the  restrictions  of  high  voltage  and  small  velocities. 
'Dmpirical  Current-Velocity  delation  for  Low  Voltages 

A  strictly  empirical  reduction  of  the  data  gives  the  simple 
representation  illustrated  in  figure  4.  This  relation  is  valid 
for  lower  applied  voltages,  l.c.,  below  the  region  of  validity 
of  the  equations  of  a  one-dimensional,  space  charge  limited  cur¬ 
rent.  The  current  is  Implicitly 


1  a  fl  +  c  ugl  ® 

i0  L  i  J 


where  c  *  constant 
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IV,  PR3S3U3B  DIFFERENCE  ACROSS  THE  DISCHARGE 
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Calculated  Treasure  Difference 

The  one-dimensional  pressure  gradient  is 


d£  a  «I3  • 
dx 


i  5 
u  +  ks 


hsing  the  simplified  expression  for  the  field  strength,  vis., 

E  a  1  [  akjx+  u?  1  ^  -  u 
v  L  J  f 


dx 

Letting  Ap  > 

Ap 


"4  [X  (ahjx+u?)^j 

*  P»«J,  “  P<«0 

"  ki  [*  (ajkx  +  u"')*]  dx 

•  il  [  1-i  up  (  !•+ ur2)^+ur2J  where  uP 


u 


VaJkT 


Eor  ur«  i 


Ap  -  li  [  1  -  i  ur] 

This  represents  the  pressure  difference  due  to  a  space 
charge  limited  current  for  airstream  velocities  small  in  re¬ 
lation  to  the  maximum  ion  drift  velocity. 

The  above  expression  for  the  pressure  involves  only  the 
current  and  the  airstream  velocity!  the  applied  voltage  does 
not  appear  explicitly.  An  expression  for  the  pressure  dif¬ 
ference  as  it  varies  with  airstream  velocity,  only;  may  be  ob¬ 
tained  by  substituting  the  exrressio-  for  the  current  as  a 
function  of  velocity  as  previously  derived. 

mhe  current  as  linearized  for  small  velocities  is 


JL 

Jo 


1+3 


aJov* 


83 


•4*- 


1 


This  gives  to  first  order  terms 


or  approximately, 

A  p  .  jUi  [  1+| 

Thie  rather  encouraging  expreeaion  indicates  that  the 
pressure  difference  increases  if  the  airstream  velocity  is 
allowed  to  increase.  The  relation  was  not  ohecked  experi¬ 
mentally  since  measurements  were  obtained  only  for  velocities 
determined  by  fixed  boundary  conditions.  As  will  be  seen  the 
liiuits  placed  on  the  velocity  by  frictional  losses  are  quite 
restrictive. 


i 
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Keasuremente  of  Induced  Velocity 

Since  the  pressure  difference  acroee  a  point  corona  die- 
charge  is  of  the  order  of  10**?  inches  of  water  at  atmospheric 
pressure,  it  was  decided  that  it  would  be  much  easier  to  mea¬ 
sure  the  induced  velocities.  To  do  this  a  hot  wire  anemometer 
wan  constructed  and  calibrated  on  a  whirling  arm. 

I.n  order  to  take  measurements  with  the  hot  wire  equipment, 
which  was  at  ground  potential,  it  was  necessary  that  the  in¬ 
duced  flow  be  transported  to  a  region  outside  of  the  influence 
of  the  discharge.  ly  extending  the  cylindrical  electrode  suf¬ 
ficiently  downstream  the  influence  of  the  hot  wire  probe  became 
negligible.  This  was  checked  by  noting  changes  in  the  corona 
current  and  by  measuring  the  discharge  current  picked  up  by  the 
probe.  Aa  a  result  of  minimizing  the  effects  of  the  probe  on  the 
discharge,  the  increased  wall  friction  is  by  no  means  small  in 
relation  to  the  pressure  forces  in  the  discharge.  Tor  example, 
the  pressure  drop  due  to  wall  friction  w -a  of  the  order  of  two- 
thirds  the  pressure  difference  across  the  discharge. 

In  figure  5  two  typical  induced  velocity  profiles  at  the 
downstream  end  of  the  tube  are  illustrated.  Included  is  a  vel¬ 
ocity  profile  without  the  discharge  showing  the  magnitude  of  the 
effect  of  the  wall  friction.  Figure  fi  is  a  plot  of  the  center- 
line  velocities  as  a  function  of  the  current.  'Vhile  there  is 
some  variability  in  the  plotted  data  the  severity  of  actual  var¬ 
iations  is  not  shown  since  only  time  average  values  of  the  vel¬ 
ocity  rre  plotted.  Velocity  differences  for  curve  (a)  in  figure 
5  v/ere  of  the  order  of  ±  .r>  ft/sec.  From  the  outset  of  the  ex- 
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pertinent  a  which  included  &  number  of  electrode  armagements 
these  velocity  fluctuation*  which  could  not  he  correlated 
with  voltage  or  current  fluctuations  were  encountered* 

Diacnseion  of  Telocity  Variations 

The  wide  ranges  of  induced  Telocity  which  had  been  attained 
for  a  fixed  applied  potential  and  constant  current  suggested 
that  the  variation  was  due  to  changes  occurring  in  the  stream  a- 
way  from  the  immediate  neighborhood  ox’  the  point.  ?or  a  fixed 
current,  conditions  at  the  point  should  he  constant.  It  was  then 
realized  that  the  time  dependency  of  the  ion  mobility  due  to 
ion-molecule  combination  might  be  of  importance.  One  way  of 
checking  this  was  to  rove  roe  the  polrity  of  the  electrodes  to 
s-.o  if  the  variations  persisted  for  the  negative  ion.  %’or  the 
atPiS  geometric- 1  arrangement  the  only  factor  which  would  <e 
different  would  be  the  mobility. 

To  accomplish  this  check  for  the  same  geometry  and  to  ko>.p 
the  hot-wire  probe  away  from  the  "hot"  electrode  the  arrangement 
shown  at  the  top  of  figure  7  was  used  with  the  hot-wire  upstream 
and  a  center  electrode  support  added.  This  support  further  add¬ 
ed  to  frictional  losses  and  gives  an  increased  amount  of  turbu¬ 
lence  at  the  point  as  compared  with  the  original  arrangement. 
There,  the  center  electrode  extended  as  a  thin  rod  from  the  stag¬ 
nation  region.  Towevcr,  these  factors  probably  do  not  detract 
from  the  primary  purpose  of  illustrating  differences  due  to  a 
difference  in  the  don  mobility. 

Tn  figure  7  the  induced  velocity  as  a  function  of  measured 
current  shows  a  range  of  velocities  for  a  given  current  for  both 


the  positive  arid  negative  point  discharge.  A  moat  important 
feature  is  that  the  minimum  induced  velocity  for  the  positive 
diacharge  is  very  nearly  the  same  as  the  maximian  for  the  neg¬ 
ative  discharge.  An  additional  point  to  note  was  that  one  com¬ 
plete  run  was  eoraplet*  d  under  approximately  the  same  conditions 
during  which  the  induced  velocity  wan  stable.  This  stable  velo¬ 
city  corresponded  to  the  minimum  for  the  positive  diacharge. 

An  additional  factor  to  be  considered  is  the  effect  on  the 
mobility  of  the  moisture  content  of  the  air.  This  not  only 
-uffects  the  value  of  the  mobility  but  the  time  dependency  as  well 
The  relative  humidity  for  the  experiment  in  question  was  meas¬ 
ured  at  47,'J  using  a  sling  paychi’ometer. 

(1) 

It  has  been  stated  that  the  mobility  of  positive  ions 
is  only  slightly  affected  by  moisture.  This  result  actually  re¬ 
fer  o  to  the  so  c.-' L led  "aged”  positive  ion.  Tyndall  and  Orind- 
(S) 

ley  found  that 

"The  initial  positive  ion  has  a  mobility  which  is  indis¬ 
tinguishable  from  that  of  a  negative  ion,  and  which  is  affected 
by  water  vapor  in  the  same  way" ■ 

"Thus,  if  the  mobility  that  is  of  importance  here  is  basically 
that  of  the  initial  positive  ion,  then  appreciable  changes  due 
to  moisture  content  can  be  expected  to  occur.  However,  they 
might  be  expected  to  be  of  Importance  only  for  the  situation 
where  the  "life11  of  the  ion  is  not  very  long  or  very  short  com¬ 
pared- to  the  characteristic  time. 
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Calculation  of  Pressure  Difference  from  Velocity  Measurements 

In  all  the  measurements  obtained  of  the  Induced  velocity 
the  results  ore  somewhat  beclouded  by  the  magnitude  of  the 
frictional  losses  and  the  variability  noted  in  the  velocity. 

For  a  best  estimate  of  the  pressure  difference  across  the  dis¬ 
charge  use  is  made  of  the  data  of  figure  7  for  the  steady  pos¬ 
itive  discharge,  i.e.  the  dividing  line  between  negative  and 
positive  discharge.  For  this  case  the  ion  mobility  appears  to 
be  "pinned  down"  at  least  to  some  constant  value.  The  pressure 
drop  due  to  frictional  losses  is  estimated  from  velocity  Burveya 
at  the  exit  of  the  tube  in  the  absence  of  the  discharge. 

From  these  velocity  surveys  the  mean  velocity  may  be  cal¬ 
culated.  "’his  is  approximately  the  velocity  at  the  entrance  to 
the  tube,  where  the  velocity  profile  Is  nearly  constant.  The 
frictional  pressure  drop  is  then  correlated  with  that  for  the 
induced  flow  at  the  same  mean  velocity.  The  pressure  drop  due 
to  friction  is  readily  obtained  by  applying  Bernoulli's  equation 
to  the  core  of  the  flow  which  is  practioally  unaffected  by  fric¬ 
tion,  i.e.,  unaffected  by  viBcous  shear,  not  unaffected  in  the 
sense  that  the  core  is  accelerated.  Thus,  the  pressure  drop  is 
obtained  using  only  a  velocity  survey  at  the  exit  of  the  tube. 

In  the  absence  of  frictional  losses  in  the  tube  the  pressure 
difference  across  the  discharge  would  be  just  equal  to  the  dy¬ 
namic  head  measured  at  the  inlet.  Actually,  the  pressure  dif¬ 
ference  due  to  the  discharge  is  the  sum  of  the  ihlet  head  plus 
the  magnitude  of  the  Treasure  drop  due  to  friction.  In  figure  8 


is  shown  the  magnitude  of  the  frictional  loss  at  the  measured 
mean  Telocity,  In  the  same  figure  the  total  pressure  differ¬ 
ence  due  to  the  discharge  1b  plotted  ss  a  function  of  the  cur¬ 
rent,  Approximately  two-thirds  of  the  preosure  increase  due  to 
the  discharge  is  loot  as  a  result  of  friction. 
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Comparison  With  Calculated  pressure  Difference 

The  data  of  figure  8  is  beet  compared  with  the  expression 
derived  for  the  variation  of  the  pressure  difference  as  a  fun¬ 
ction  of  the  current  and  the  velocity  ratio*  Ur,  vis* 

Ap  »  £l  -  £  ur] 

The  limitations  imposed  on  the  velocity  which  are  primarily 
due  to  frictional  losses  give  a  value  of  U?  of  the  order  of 
0*02*  Thus*  the  effect  of  the  airstream  velooity  gives  cor¬ 
rections  of  the  order  of  1$  which  may  be  neglected* 

If  the  axial  distance  between  the  electrodes  is  taken  as 
l  and  the  current  density  computed  using  the  cross-sectional 

■*> 

■.*.  area  of  the  tube,  then  the  value  of  k  which  matches  the  straight 

line  of  figure  8  is  k  a  ?*4  x  10“3  ft/sec  •  This  is*  to  two 

vo  It  is/ ft 

significant  figures,  the  value  of  k  given  on  page  6  and  correc¬ 
ted  for  the  air  density  of  the  experiment.  It  is  probably  the 
correct  value  to  use  since  the  "life"  of  the  ions  is  of  the 
order  of  10“ 3  sec.,  and,  especially  with  moisture  present,  the 
mobility  should  correspond  to  the  value  of  that  for  "new"  pos¬ 
itive  ions* 


y._  .coronjsioHS 


1. )  The  restrictions  placed  on  the  induced  velocity  preclude 
confirmation  of  the  equations  for  the  pressure  difference  to 
anywhere  near  the  limits  for  whioh  the  characteristics  of  the 
discharge  in  an  airstream  were  obtained.  However*  the  data 
indicates  that  the  one-dimensional  analysis  gives  essentially 
the  correct  result  for  low  speeds  for  the  pressure  difference 
and*  in  addition*  gives  a  reasonable  description  of  the  cur- 
rent-velocity  relation  for  fairly  high  speeds.  Thus,  predic¬ 
tions  based  upon  this  model  should  not  be  too  far  from  the 
truth. 

2. )  This  study  was  intended  to  raiee  questions  regarding  ap¬ 
plications  of  the  basic  process.  The  relation  of  the  ion  mo¬ 
bility  with  the  properties  of  atmospheric  air  and  the  geometry 
of  a  particular  problem  is  juet  such  a  question  which  presented 
itself  and  which  was  not  answered. 

5.)  With  respect  to  aeronautical  engineering  applications  of 
electrical  pressure  forces,  it  waB  realised  from  the  outset  that 
this  study  might  be  classified  in  the  realm  of  "bird  -  watching”. 
Power  inputs  are  email  and  power  outputs  in  terms  of  directed 
energy  are  etlll  much  smaller. 

The  efficiency  of  the  experimental  device  ae  a  pump  can  be 
estimated  from  the  expressions  for  the  pressure  difference  and, 
for  a  given  applied  voltage,  the  expression  for.  the  current  ae 
a  function  of  the  airstream  velocity.  The  power  input  is  Juet 
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Vi  and  the  power  output  is  uAAp.  As  an  example,  from  the 
data  of  figure  1  for  a  fixed  applied  voltage  of  11,600  volts, 
i0  ■  J?«r2/*amps,  and  considering  airstream  velocities  email  in 
comparison  to  the  maximum  ion  drift  velocity  the  efficiency  in 
percent  is 

m  ,P?u  %  (u  in  ft/sec} 

For  a  mean  induced  velocity  of  the  order  of  5  ft/sec.  the  ef¬ 
ficiency  of  the  device  as  a  pump  is  of  the  order  of  1  $£, 

The  efficiency  is  limited  to  low  values  for  two  reasons* 
a.)  The  tpobi:  important  restriction  is  that  placed  on  the 
velocity  of  the  flow.  This  is  primarily  the  result  of 
frictional  losses. 

h.)  The  corora  discharge  is  an  inefficient  source  of  ion¬ 
ization,  i.e.,  the  useable  power  is  (V-VB)i  whereas  the 
total  power  input  is  Vi.  V8  is  of  the  order  of  thou¬ 
sands  of  volts  whereas  the  useful  ionisation  requires  an 
ionizing  potential  only  of  the  order  of  15  volts. 

To  increase  the  magnitude  of  the  effects  considered  here, 
J.e.,  to  introduce  much  higher  charge  and  current  densities 
other  methods  are  certainly  available.  For  example,  the  mech¬ 
anism  of  ionisation  might  employ  alpha  rays  rather  than 
the  corona  discharge.  The  magnitude  of  the  effects  due  to  a  cor¬ 
ona  discharge  ls  determined  by  the  arcing  limit.  The  high 
field  intensity  required  for  the  sustenance  of  the  discharge 
limits  the  current  hclow  the  disruptive  discharge  to  relatively 
low  values. 
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